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We investigate the effect of hyperfine interactions on ultracold molecular coUisions in magnetic fields, 
using ^^Mg(^S) + ^''NH(^S~) as a prototype system. We explore the energy and magnetic field 
dependence of the cross sections, comparing the results with previous calculations that neglected 
hyperfine interactions (Phys. Rev. Lett. 103, 183201 (2009)). The main effect of hyperfine interac- 
tions for spin relaxation cross sections is that the kinetic energy release of the dominant outgoing 
channels does not reduce to zero at low fields. This results in reduced centrifugal suppression on 
the cross sections and increased inelastic cross sections at low energy and low field. We also analyze 
state-to-state cross sections, for various initial states, and show that hyperfine interactions introduce 
additional mechanisms for spin relaxation. In particular, there are hyperfine-mediated collisions to 
outgoing channels that are not centrifugally suppressed. However, for Mg-|-NH these unsuppressed 
channels make only small contributions to the total cross sections. We consider the implications of 
our results for sympathetic cooling of NH by Mg and conclude that the ratio of elastic to inelastic 
cross sections remains high enough for sympathetic cooling to proceed. 

PACS numbers: 34.50.Cx, 33.15.Pw, 37.10.Mn, 37.10. Pq 



I. INTRODUCTION 

High-density samples of cold (T < 1 K) and ultra- 
cold (r < 1 /iK) molecules are likely to have impor- 
tant applications in fields including quantum informa- 
tion science, high-precision spectroscopy and quantum- 
controlled chemistry [I]. For nearly three decades, ul- 
tracold atoms have been produced using laser Doppler 
cooling and evaporative cooling. At temperatures be- 
tween 1 nK and 1 /iK, atoms enter a fully quantal regime 
with novel properties, forming Bose- Einstein condensates 
[2H1] and Fermi-degenerate gases E] ■ However, cooling 
molecular samples to the temperatures required for quan- 
tum degeneracy is much more difficult than for atoms, 
because molecules have a much richer internal structure. 

Molecular cooling techniques can be classified as ei- 
ther direct or indirect. In direct methods, molecules are 
cooled from relatively high temperatures by techniques 
such as buffer-gas cooling fTl, Stark deceleration 8] or 
Zeeman deceleration [9 llj. In addition, laser Doppler 
cooling has very recently been demonstrated for SrF [H] , 
although only a very limited number of molecular species 
are likely to be amenable to this technique. Indirect 
methods, by contrast, form ultracold molecules in gases 
of previously cooled atoms by magnetoassociation [OlfTi] 
or photoassociation [13l[T5]. In the last few years, indi- 
rect methods in association with stimulated Raman adi- 
abatic passage (STIRAP) |Jjj have succeeded in produc- 
ing ground-state molecules at temperatures below 1 
[17l - fT9| . Indirect methods are however restricted to the 
relatively few atomic species that can be laser-cooled. 

Direct methods can in principle be used for a wide va- 
riety of molecules, but the lowest temperatures achieved 
so far are in the millikelvin regime. A second-stage cool- 
ing technique is needed to cool molecules down to the 
microkelvin range, from which evaporative cooling could 



reach the region of quantum degeneracy. Sympathetic 
cooling is a promising second-stage cooling method, 
which relies on thermalization of the molecular species 
of interest by collisions with ultracold atoms. Sympa- 
thetic cooling has been used for almost three decades 
to cool molecular ions to temperatures around 100 mK 
and below [20, ch. 18], and has more recently also been 
used for neutral atoms [211 ES] . Its extension to neutral 
molecules was first proposed in 2004 by Soldan and Hut- 
son [221 . Since then, accurate quantum scattering calcu- 
lations have suggested that it should work for a few di- 
atomic [231 ES] and polyatomic molecular species [26l [27] 
in combination with appropriate ultracold atoms. Al- 
though experimental validation is still needed, there is 
great hope that sympathetic cooling will eventually be- 
come a 'routine' second-stage cooling technique, capable 
of producing large molecular samples in the microkelvin 
regime, where most of the new interesting applications 
would become possible [1]. 

Cold atoms and molecules are usually held in traps 
formed with electric, magnetic or optical fields. By far 
the most experimentally accessible traps are those formed 
by static electric [55] and magnetic fields [7]. The main 
limitation of such traps is that they can hold only those 
molecules whose energy increases with applied field, i.e. 
those in low-field-seeking states. However, there are al- 
ways untrapped high-field-seeking states that lie energet- 
ically below the low-field-seeking states. Inelastic (deex- 
citation) collisions that transfer molecules to high-field 
seeking states therefore lead to trap loss. Elastic colli- 
sions, by contrast, are required for thermalization. Thus, 
the success of both sympathetic and evaporative cooling 
depends on a favorable ratio of elastic to inelastic cross 
sections, of at least about 100 [^ . 

In the (ultra)cold regime, hyperfine interactions are 
larger than or comparable to the kinetic energies in- 
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volved. Hyperfine couplings also provide both additional 
relaxation mechanisms and new ways to control atomic 
and molecular collisions. However, only a few calcula- 
tions have so far considered the consequences of hyper- 
fine interactions on cold molecular collisions. Bohn and 
coworkers considered the scattering of polar ^^OH(^n3/2) 
molecules in static electric 30 and magnetic fields [3T| . 
Lara et al. [32} ,33j studied the field-free scattering of 
Rb(2S) and OH(2n3/2). Tscherbul et al. studied (ul- 
tra)cold collisions of YbF(^E) molecules with He in ex- 
ternal electric and magnetic fields [Mj. They found 
that simultaneous electron and nuclear spin relaxation 
for collisions in the ground rotational state can occur 
through couplings via excited rotational states and hy- 
perfine terms. 

In the present paper, we explore the effects of hyper- 
fine interactions on the ultracold scattering of NH with 
Mg in the presence of a magnetic field. This extends the 
study in [M] in various ways: (1) NH is a instead of a 
molecule; (2) NH has two nuclei with non-zero spin, 
leading to more complicated hyperfine structure than for 
YbF; (3) we consider hyperfine states other than spin- 
stretched states. Our results are compared with those of 
a previous study of Mg-NH neglecting hyperfine terms 
[24j . and allow us to make some general conclusions on 
how hyperfine interactions can affect the prospects of 
sympathetic cooling. 

II. THEORY 

We consider the scattering of ^'*NH(^S~) molecules 
with 24Mg(iS) atoms in the presence of an external 
magnetic field B, whose direction defines the space- 
fixed Z-axis. The system is described using Jacobi vec- 
tors/coordinates: the vector r runs from N to H while 
R runs from the center of mass of NH to Mg. The angle 
between r and R is 9. The collision is studied by solving 
the time-independent Schrodinger equation for the scat- 
tering wave function at energy E, 7i'i> = E'i', where 
H is the Hamiltonian for the colliding pair. 

A. Effective Hamiltonian 

By convention, lower-case symbols are used for quan- 
tum numbers of the individual monomers, and capital 
letters are used for quantum numbers of the collision 
complex as a whole. Where necessary, the subscripts 
1 and 2 refer to Mg and NH respectively. The diatom is 
considered to be a rigid rotor in its ground vibrational 
state. The effective Hamiltonian can be written 

+ (1) 

where L is the space-fixed operator for the end-over-end 
rotation and fi is the reduced mass of the complex. In 



general, T^mon contains terms describing both isolated 
monomers, 'Hmon = "Hi -t- 'H2- However, ground-state 
^^Mg has both zero electron spin and zero nuclear spin, 
so that it contributes only a constant energy and it is 
convenient to set Hi — 0. Finally, Hu includes all in- 
teractions between the monomers, which in the present 
case reduces to the potential energy surface V{Ur,R), 
conveniently expanded in Legendre polynomials 

V{ur,R) = £ Vk{r,^,R)Pk{cose). (2) 

k=0 

Here Ur denotes a unit vector in the direction of r, r^q is 
the equilibrium distance of NH in its ground vibrational 
state and are the radial strength functions. The po- 
tential for the Mg-NH system was reported in Ref. ^35, . 
The Hamiltonian for an isolated NH('^I]~) molecule 

can be written "^2 = ?irot+?isn+^ss+^hf +?iz- The dif- 
ferent terms will be discussed below and correspond, re- 
spectively, to the rotational, electron spin-rotation, elec- 
tron spin-spin, hyperfine and Zeeman interactions. 

If centrifugal distortion and all other higher-order cor- 
rections are neglected, the Hamiltonian for the mechan- 
ical rotation of NH is Hrot = ^nh^^, where 6nh is the 
rotational constant (with dimensions of energy) and n is 
the operator for the rotational angular momentum. The 
electron spin-rotation term, arising from the interaction 
between the magnetic moment associated with the com- 
posite electronic spin of NH, s, and the magnetic field 
induced by its rotation, can be written Hsn = in ■ s, 
where 7 is the spin-rotation constant. 

The direct dipolar interaction between the unpaired 
electrons in NH('^I]~) may be written [36 

iiss ~ 2AssT^^o(J, s)T^^o(itr, Ur), (3) 

where Ass is the electron spin-spin constant and T*"' rep- 
resents a spherical tensor of rank k, with g-component 

<?■ 

Since both the ^^N and ^H nuclei have non-zero nuclear 
spin, jn = 1 and in = 1/2, the hyperfine Hamiltonian can 
be written 

Hm — 'Hsi.N+'HF.N+'Hin,N+'HQ,N+'Hsi.H+'HF.H+'Hin,H. 

(4) 

Here, Hsi represents the direct dipolar interaction be- 
tween the magnetic moments associated with a given nu- 
clear spin I and the composite electron spin of the open- 
shell S molecule and can be written 

?isi = -^/T0g5/^Bg^m(A*o/4^)Tl(s,C2) •Ti(J) 

« V6toT2^o(s,«), (5) 

where gs and gi are the electron and nuclear ^-factors, 
/iB and /iN are the Bohr and nuclear magnetons, and 
/io is the magnetic permeability of free space. The axial 
dipolar interaction parameter is related to the widely- 
used constant c defined by Frosch and Foley [37] as — 
c/3. 
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The Fermi (or Ferrai-Breit) contact interaction Hp oc- 
curs whenever there is a non-zero electron-spin density 
at a nucleus with non-zero spin i. It may be written 
Jip = bps ■ i, where bp is a coupling constant given 
by bp = (2/3)<7s^B5«mlV'(0)P, with |V'(0)|2 the spin 
density. The coupling constant can also be written 
in terms of Frosch and Foley's b and c parameters as 
bp = b+{l/3)c. 

Since > 1/2, the interaction between the nuclear 
electric quadrupole moment Q and the electric field gra- 
dient at the ^^N nucleus V-E must be included. In gen- 
eral, this is 



nQ = -eT^{VE)-T^{Q), 
and for a diatomic molecule reduces to 



eqoC 



4i{2i - 1) 



(6) 



(7) 



where qo is the electric field gradient. 

The nuclear spin-rotation terms 'Hin are the nuclear 
counterpart of the Hgn interaction discussed above, and 
are written "Hin — Cifi ■ i, where Ci is the corresponding 
nuclear spin-rotation constant. 

Finally, if only electron and nuclear spin Zeeman terms 
are taken into account, 

"Hz = ffsMBSzJB - /iN ^ 5x«XzB(l - o-x). (8) 

X=N, H 

The nuclear shielding factors cx are extremely small and 
are neglected in the present work. 

The molecular constants used for the -'^"'NH('^E^) rad- 
ical are listed in Table HI 



B. Coupled-channel equations 

We solve the scattering problem using the coupled- 
channel method. First, the total wave function is ex- 
panded in a set of N conveniently chosen basis functions 



(9) 



Here, f is a collective variable including all coordinates 
except R, and a is the set of quantum numbers that la- 
bel our basis functions. Each different combination of 
quantum numbers a is said to define a channel. A set of 
coupled differential equations for the channel Junctions 
Xa{R) is then obtained by substituting 'I'(i?,^) into the 
time-independent Schrodinger equation. 



Xa 



dR^ 



E 



{Waa' - e.5aa')Xa' 



(10) 



where 5ij is the Kronecker delta, e = IpuEjf? is a scaled 
energy, and 



Waa'{R) 



2^ 



2^i?2 



(11) 



The coupled equations ( 10 1 are solved by propagating a 



complete set of independent solution vectors from i?niin) 
deep in the inner classically forbidden region, to i?max, 
large enough that the effects of the interaction potential 
have died off. If necessary, the solutions are transformed 
at i?max into a basis set in which Waa' and are diago- 
nal at i? = 00 131] , and the transformed channel functions 
are matched to the standard scattering boundary condi- 
tions [42]. This gives the scattering matrix S, from which 
all quantities of interest, such as state-to-state cross sec- 
tions and scattering lengths, may be calculated. Numer- 



ical details are given in Sec. IIIB 



Basis set and matrix elements 



TABLE I: Molecular parameters for "NH(^E", v = 0). 



Parameter Value Reference 

feNH/cm-^ 16.343 |31j 

7/cm-^ -0.055 [38] 

Ass/cm-i 0.92 [58] 

gN 0.40376 [39] 

&f,n/MHz 18.83 @D] 

cn(= 3to,N)/MHz -67.922 gO] 

(egoQ)N/MHz -2.883 [iQ] 

C7i,n/MHz 0.1455 gD] 

gn 5.58568 [39] 

6f,h/MHz -66.131 [iO] 

ch(= 3to,H)/MHz 90.291 [ID] 

C7i,h/MHz -0.061 [40] 



We use a fully uncoupled basis set \a) = \a) \LMl), 
where \a) = \i-NiTiiN) I^HWin) \sms) |nm„) describes the 
state of the monomers, and niA (or Ma) denotes the 
projection on the field axis of the vector operator A. 
None of the interactions considered here change the elec- 
tronic or nuclear spins, so we omit the labels s and i and 
label our basis functions {mM,miii,ms,n,mn, L, AI^). 
A static magnetic field conserves both the projection 
Mtot of the total angular momentum F and the total 
parity P of the system. These are explicitly Mtot = 
TOiN + "T-iH + "T-s + w„ -t- Ml and P = pip2(— 1)^, with 
Pi — I the parity of Mg(^S) and P2 — (—1)"+^ the parity 
of NH('^I]^). The matrix elements for the centrifugal, ro- 
tational, electron spin-spin and interaction potential are 
diagonal in, and independent of, the nuclear spin quan- 
tum numbers. The corresponding expressions in our basis 
set are thus readily obtained from Ref. [¥T| . 
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The Zeeman matrix elements are completely diagonal in the uncoupled basis set, 

(sTOsI (jhTOihI («n"i,:n| 'Hz |iN™,:N) |iH"^ffl) \s'm^) = B [gsfJ^sms - (ffN'TiiN + 5hto^h)] • (12) 

Here and throughout this section, the matrix elements are fully diagonal with respect to quantum numbers that do 
not appear explicitly in the expression. 

The Fermi contact term and the electron and nuclear spin-rotation terms all share a similar structure, "Hjj^jj = kJi ■j2, 
where k is a scalar while ji and j2 are vector operators. In general, their matrix elements in a decoupled basis set, 
\jimji) \j2mj2), are 

+'5m,im;.,±l^m,2™^2Tl2 + [^2 O2 + 1) - mjsm^a] • (13) 

Such terms can mix functions with adjacent values of the projections of ji and j2j but preserve the sum toi2 = 
niji + mj2- 

The electron- nuclear spin dipolar interaction, Eq. ([5]), has matrix elements 

(nm„| {srusl {im^lHsi |"'"4) 

= io\/30(-l)''"'"'+'""'=^"" + l)(2i + l)s{s + l)(2s + l)(2n + l)(2n' + 1)]^^^ 

^fn2n'\y^fll '2 \ f i 1 *W * ^ " ^ ^' \ (14) 

where q = qi + q2 and (:::) is a 3-j symbol. This term produces couplings off-diagonal in one nuclear spin projection 
rrii, along with rus and m„ (keeping their sum unchanged), with An = n' — n = or ±2. This latter selection rule is 
required to conserve p2 ■ 

Finally, the quadrupole interaction for the ^^N nucleus has matrix elements 

/ I/- l-i) I- ' \ l ' '\ (e9o'9)N. 1^^^,-mi^^-m„ r/o , ' , imVz/ U 2 v! \ ( In 2 In \ 

(nm„|(zNm,N|^Q,N|«N"i,N) = (-1) " [(2n l)(2n + 1)] ' n o n • o ■ 

4 \ / \ —In In / 



X 

p 



Yi-lf i n 2 n' \l ^n 2 *n ^^^^ 



These couple functions An = 0, ±2 and Ato^n = — Ato„ = 0,±1,±2, thus preserving the sum to^n + "nin as well as 

1 



The first 3-j symbol in Eqs. ( 14 ) and ( 15 ) implies that 



the electron-nuclear spin and quadrupolar interactions 
have no direct off-diagonal matrix elements between n = 
functions, so that their dominant contribution for n — 
is a second-order effect necessarily involving the n = 
2 excited rotational state; this is also the case for the 
electron spin-spin interaction. 



III. RESULTS AND DISCUSSION 

A. Hyperfine levels 

The energy levels correlating with the ground rota- 
tional state of an isolated ^^NH('^E^) molecule, in an 
external magnetic field, are shown in Fig. [l] The solid 
lines correspond to the inclusion of all hyperfine inter- 
actions as described in Sec. Ill A| while the dashed lines 



show the energy levels obtained when hyperfine terms are 
neglected. 

The hyperfine-free levels are labeled at zero field by 
the eigenvalues of the angular momentum j — n -f s, 
which has only one allowed value, j = 1, for n = and 
s = 1. In the presence of a magnetic field, a level with 
quantum number j splits into 2j -\- 1 Zeeman components 
characterized by the projection rrij onto the field axis, 
and the corresponding eigenstates can be represented as 
l(ns)jmj). For n 
alternatively be labeled by 

The pattern is much more complicated when hyper- 
fine interactions are included. In this case, the zero-field 
levels are labeled by the eigenvalues of the total angular 
momentum / resulting from coupling the nuclear spins 
«N and iH with fi and s. In general, three or more angu- 
lar momenta can be coupled using a variety of schemes. 
We start, as before, by first coupling fi and s to form 
J. Then, given the molecular constants in Table [ij it is 



0, rrij — m., and these states can 
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FIG. 1: (Color online) Magnetic field dependence of the en- 
ergy levels of ^''NH('^E~) correlating with the ground ro- 
tational state, including (black, solid) and excluding (blue, 
dashed) hyperfine terms. 



convenient to couple j first with in to produce a resul- 
tant /fj. Finally, /h is coupled with in to give /. In 
the particular case of n = 0, j = s = 1 and /h = 1/2, 
3/2, which produces levels with / = 1/2 and 3/2 for 
/h = 1/2 and / = 1/2, 3/2 and 5/2 for /h = 3/2, as 
shown in Fig. [T] For rt = 0, the Fermi contact terms 
make by far the largest contributions to the splittings, 
and levels with /h = 3/2 lies below those with /h = 1/2 
because 6f,h < 0. In a magnetic field, each / state fur- 
ther splits into 2/ -t- 1 sublevels, producing a total of 
(2iN + l)(2iH + l)(2s -t- 1) = 18 components correlating 
with the ground rotational state. At low fields, below 
about 10 G, the eigenfunctions are approximately repre- 
sented as |(ns)i, (i«H)/H, (/h*n)/"^/)- In what follows, 
even though /h and / are not good quantum numbers in 
a field, we use quantum numbers (/h, /, ^ti/) in parenthe- 
ses to identify the states at low field. In addition, levels 
will be labeled fii (i = 1, 18) in order of increasing energy 
at fields above 50 G, where TOs is a nearly good quan- 
tum number. With this convention, correspond to 
nis = -1, /37-/3i2 to rris = 0, and /3i3-/3i8 to = +1. 

For isolated NH, the total projection is a good 
quantum number. Hence, as a function of field, states 
corresponding to different mf can cross while states of 
the same mf cannot. However, in our model Hamilto- 
nian, states /3io and /3i3 (both with — 1/2) are seen to 
cross at about 25 G. This is a nonphysical effect which re- 
sults from neglecting the interaction between the nuclear 
spins of N and H, usually written in the form C4£n • in- 
For NH, C4 is extremely small and has not been measured 
experimentally; there is in reality an avoided crossing be- 
tween states /3io and /3i3, but it is extremely tight. 

At high field, the terms in the monomer Hamiltonian 
that mix states with different values of nis are small with 
respect to the electron Zeeman splitting and nis is well 
defined. For n = at fields over 75 G, three groups of 
NH levels corresponding to rn^ = —1,0 and -1-1 can be 



identified, containing six hyperfine levels each. In this 
regime, the individual nuclear spin projections uim and 
m^H are also nearly conserved and the eigenfunctions are 
well represented by individual basis functions \a). Above 
75 G, more than 93.5% of any eigenstate is represented 
by a single basis function \a). In the high-field limit we 
label states with quantum numbers in square brackets, 

[ms,TOiN,miH]. 

All TTLs = +1 states are low- field-seeking and there- 
fore trappable in a static magnetic trap. When hyperfine 
interactions are included, there are six such levels with 
different values oi mf. In contrast, when hyperfine terms 
are neglected, there is only one such state. 

One state of particular interest is the spin- stretched 
state, in which / and mf take their highest possible 
values. Except for terms off-diagonal in n, this state 
is exactly represented by a single basis function in ei- 
ther possible basis set, (/h,/, "i/) = (3/2, 5/2, +5/2) or 
[ms,TOiN,miH] = [+l,-)-l,+l/2]. For "NH, the spin- 
stretched state is /3i5, and lies below three levels from 
the (/h,/) = (1/2,3/2) manifold, /3i6~/?i8. 



B. Scattering cross sections 

We have carried out scattering calculations using the 
MOLSCAT package [13], as modified to handle colli- 
sions in external fields [H]. The coupled equations were 
propagated with the hybrid log-derivative Airy method 
of Alexander and Manolopoulos using a fixed-step- 
size log-derivative propagator for 2.5 < i? < 50 A, with 
AR = 0.025 A, and a variable-step-size Airy propaga- 
tor for 50 < i? < 250 A. To a good approximation, the 
computer time is dominated by operations on relatively 
large matrices that scale with the total number of chan- 
nels TV as c!c N^. The time needed to perform a cal- 
culation including hyperfine interactions is thus approxi- 
mately [(2iN + l)(2iH + 1)]^ = 216 times larger than that 
required for an equivalent calculation neglecting them. 
In order to make our calculations tractable, the basis set 
used in the present work was reduced slightly from that 
used in Ref. [23], to rimax — 5 and Lmax = 6. Under 
these conditions iV « 1, 500, with the actual number de- 
pending on the initial state, Mtot and P. The reduction 
in n-max and imax docs not change the results by more 
than about 5%. 

NH molecules that undergo transitions between hyper- 
fine levels of the TOs = +1 manifold remain in a magnet- 
ically trappable state. However, the associated kinetic 
energy release ranges from 0.7 to 5.8 mK and will ei- 
ther heat the trapped gas or eject one or both collision 
partners from the trap. We will thus assume that all in- 
elastic processes have a negative impact on the success 
of sympathetic cooling. In any case, as will be shown in 
Sec. jlHB 3 a[ transitions with Am.; = do not contribute 
appreciably to inelasticity. 
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1. General considerations 

a. Analytical model. The total inelastic cross sec- 
tion may be decomposed into partial-wave contributions, 



C/3,inol 



J2 ^^i- 

I3'^I3.L,L' 



(16) 



When inelastic scattering is weak compared to elastic 
scattering, the first-order distorted-wave Born approxi- 
mation [45^ provides relatively simple expressions for the 
off-diagonal 5"- matrix elements. Volpi and Bohn [35] gave 
an analytical formula for the threshold behavior of the 
partial inelastic cross sections under these conditions. 



(17) 

Here, cr^^, is a factor independent of the collision energy 
E, while AEppi is the kinetic energy released in the tran- 
sition from the state /3 to f3' . This formula was used in 
interpreting the energy and magnetic field dependence of 
hyperfine-free scattering cross sections in Mg -I- NH 

When hyperfine terms are neglected, I^Eppi in Eq. ( 17 ) 
is simply —gsUBBAnis and a^^, is independent of the 
magnetic field strength B. However, when hyperfine 
terms are included, the B-dependence of AEpi^i is much 
more intricate, as seen in Fig. [l] In particular, for many 
transitions, AEp/^i does not approach zero at low mag- 
netic fields. This has substantial effects below about 
50 G. At higher fields, AEp^r approaches its hyperfine- 
free form even when hyperfine interactions are included. 
In addition, a^^, varies with B because the character of 
the NH eigenstates depends on the magnetic field. 

b. Relaxation mechanisms. The quantum- 
mechanical theory of electron spin relaxation in 
collisions of '^S molecules with structureless atoms, 
neglecting hyperfine effects, was developed by Krems 
and Dalgarno [47\ In general, inelastic collisions are 
driven by the anisotropy of the interaction potential, but 
this does not have matrix elements that are off-diagonal 
in mg. Collisions that change thus occur only 
because it is not strictly a good quantum number, and 
basis functions with different values of rrig are mixed 
by terms in the monomer Hamiltonian. In the absence 
of hyperfine interactions, the only such terms are the 
spin-spin and spin-rotation Hamiltonians. 

In the NH('^S]~) case that we study here, Ass ^ 7 and 
the electron spin-spin terms are dominant. There are no 
matrix elements of Uss between n = states, because 
nin cannot change from 0. In the absence of hyperfine 
interactions, the effect of Has is to mix into the ground 
state {ms,n,mn, Ml) = (tos, 0, 0, M^) a small amount 
of rotationally excited functions (m^ — q,2,q, Al^) with 
q — 0,±1,±2. The potential anisotropy (principally V2) 
can then drive transitions between these mixed states. 
Spin relaxation for ri = states thus proceeds mainly 
via the combination of the spin-spin interaction and the 
potential anisotropy, and leads to transitions with Am^ 



compensated by a change in Ml to conserve Mtot- If 
L = 0, this requires L' > and the corresponding cross 
sections are therefore centrifugally suppressed at low en- 
ergies by barriers in the outgoing channels. In what fol- 
lows, we will refer to this as the 771am relaxation mecha- 
nism. 

When hyperfine coupling is included, states with dif- 
ferent values of and are strongly mixed at low 
fields. For the n = states, this mixing is almost en- 
tirely due to the Fermi contact interactions. Under these 
circumstances the main mechanism can drive all possible 
transitions among the n — states, although with vary- 
ing degrees of centrifugal suppression as described be- 
low. However at high fields, once the Zeeman splittings 
are large compared to the Fermi contact interactions, rus 
and rui become nearly good quantum numbers. Since the 
main mechanism does not affect the nuclear spin projec- 
tions, there is a propensity for transitions with Ami — 0. 



2. Dependence on collision energy 

Fig. [2] shows the elastic and total inelastic cross sec- 
tions as a function of energy at magnetic fields _B = 1, 10 
and 100 G, for collisions starting in the spin-stretched 
state, /3i5 ^ {fHj,mf) = (3/2, 5/2, -f 5/2), and the 
two highest-lying states, = (1/2,3/2,-1-1/2) and 

/3i8 = (1/2, 3/2, -t-3/2). At very low energies, the scatter- 
ing of incoming partial waves with L 7^ is suppressed 
by centrifugal barriers; our calculations include all con- 
tributions from incoming s, p and d waves (L = to 2), 
which gives integral cross sections that are converged at 
collision energies up to about 40 mK. The dashed lines 
in Fig. [2] show pure s-wave cross sections, and it may be 
seen that collisions with L > become significant above 
E w 10-" K. 

The dependence on the magnetic field will be analyzed 
in detail in the next section. However, in general terms 
it is clear that at high fields (100 G and above), the cross 
sections including hyperfine interactions are quite similar 
to those from hyperfine-free calculations, while at lower 
fields they are very different. In particular, the suppres- 
sion of inelastic cross sections that occurs in hyperfine- 
free calculations at low fields and low energies is much 
reduced (resulting in larger inelastic cross sections) when 
hyperfine interactions are included. 

The major effect of hyperfine interactions is that, at 
low fields, they increase the kinetic energy release and 
thus reduce the centrifugal suppression of the inelastic 
scattering. At 1 G, for example, the /3n and /3is states 
can relax to channels with kinetic energy releases of up 
to 7.2 mK, while the corresponding value for the spin- 
stretched /3i5 state is only 2.6 mK. When hyperfine cou- 
pling is excluded, however, the kinetic energy release at 
1 G is only about 270 /iK. 

For simplicity, let us consider the case of s-wave scat- 
tering. The projection of the total angular momentum, 
Mtot = iTT-f + Ml (with ruf = rris + m^N + rnm) is con- 
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10' 




Collision energy, E/K 



FIG. 2: (Color online) Elastic and total inelastic cross sec- 
tions as a function of collision energy, for various magnetic 
fields and initial states. Solid lines include s, p and d-wave 
contributions, while dashed lines are s-wave cross sections. 
Dotted lines show power-law to guide the eye. 



served in a collision, so L' > |M£| = |Ato/|. In addition, 
conservation of parity requires L' to be even. For scat- 
tering from the spin-stretched state, /3i5, the main relax- 
ation channels have = 1/2 and 3/2, for which L' must 
be at least 2. Channels with L' = 2 dominate at low ener- 
gies, because d-wave centrifugal barriers (height 23 mK) 
are much lower than (7- wave barriers (height 140 mK). 
Transitions to levels for which L' = 2 is not possible. 



such as (m'pL'^iJ = (-1/2,4), (-3/2,4), (-5/2,6) have 
negligible contributions at the magnetic fields considered 
here. 

For molecules that are initially in a non-spin-stretched 
state, inelastic collisions with Aruf = are possible. 
There are then relaxation channels with L' — 0, which 
are not centrifugally suppressed. However, transitions to 
these states are made possible only by invoking hyperfine 
couplings. Such transitions make very little contribution 
to the total inelastic cross section, which remains domi- 
nated by the (centrifugally suppressed) main mechanism. 
Even for molecules in non-spin-stretched states, the main 
effect of hyperfine interactions is through an increased ki- 
netic energy release that helps overcome the centrifugal 
barriers at low fields and low energies. 

In general, the elastic scattering depends on the phases 
of diagonal elements of the S'-matrix, which are only 
very slightly affected by the inclusion of hyperfine terms. 
Therefore, as seen in Fig. [2] the elastic cross sections 
including hyperfine interactions are very similar to the 
hyperfine-free results at all energies and fields. 



3. Dependence on magnetic field 

a. State-to-state s-wave cross sections. Fig.lSlshows 
the state-to-state s-wave inelastic cross sections o^Z^^/ = 
aois^L'p' as a function of magnetic field B, for col- 
lisions starting in the spin-stretched state, /3i5, with 
quantum numbers (/ui/j'Tiy) = (3/2,5/2,-1-5/2) and 
[jTis, TOiN, JTiiH] — [+1, +1, +1/2]. There are six main 
contributions, all to channels with L' — 2, corresponding 
to /3'4 (Am, = -2), /?^, and /3[o (Am, = -1), and 
and /3^4 (Am, = 0). The corresponding final-state 
energies are shown color-coded in the bottom panel. 

At low fields, where the main mechanism can drive 
all possible transitions, the state-to-state cross sections 
are governed by the kinetic energy release. The largest 
cross section at E ^ 10"^ K is to = (3/2, 1/2,-^1/2), 
closely followed by and which are (3/2, 3/2, +3/2) 
and (3/2,3/2,-1-1/2) respectively. These channels have 
the largest kinetic energy release and therefore experi- 
ence less centrifugal suppression. The relatively minor 
channels /3[^ and P'^ have zero kinetic energy release at 
low field, while /3[q is energetically accessible only at fields 
above about 15 G (slightly dependent on E), as seen in 
the bottom panel. 

At high fields, where m, and rrii become nearly good 
quantum numbers, transitions with Am^ = are favored. 
The two strongest channels are fi'^ = [—1,-1-1,-1-1/2] 
(-90%) and = [0,-^^+1/2] (-10%). The former 
is stronger because of the larger kinetic energy release 
associated with Am, ~ —2. The largest cross sec- 
tions to channels with Am^ ^ are those to P[q = 
[0, -hi, -1/2] and = [0, 0, +1/2], with the former mak- 
ing a greater contribution because the Fermi contact in- 
teraction is stronger for H than for N. Transitions to 
/J^g = [-^1,-1,-1-1/2] and /3[4 = [-hi, 0, +1/2] are weak 
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= (1/2, 1/2, +1/2) 
= (3/2, 1/2, +1/2) 
= (3/2, 3/2, +3/2) 
Pi, = (3/2, 3/2, +1/2) 
' = (3/2, 5/2, +3/2) 
= (3/2, 5/2, +1/2) 



,L.O 
p, inel 




Magnetic field, B/G 

FIG. 3: (Color online). State-to-state s-wave inelastic cross 
sections for collisions originating in the spin-stretched state 
/3i5, as a function of magnetic field, for collision energies of 
10"'' K (top panel) and 10~® K (center panel). The bottom 
panel shows the initial and final states energies, color-coded as 
for the cross sections. Solid, dashed and dotted lines represent 
states with rris = —1,0 and +1 respectively. 



both because of the change in to^^n and because they 
do not change and thus have a small kinetic energy 
release. 

The relative state-to-state cross sections are fairly in- 
sensitive to the collision energy, as seen by comparing 
the top and middle panels in Fig. |3] The only quali- 
tative difference is in the cross sections to /JJg and f3[^, 
for which the kinetic energy release is zero at zero field. 
The inelasticity to these states is the most affected by 
centrifugal suppression. The outgoing kinetic energy for 
these channels is thus dominated by E, and the increase 



= (1/2, 3/2, +1/2) 
= (1/2, 3/2, -1/2) 
= (1/2, 1/2, -1/2) 
= (1/2, 1/2, +1/2) 
= (3/2, 5/2, +5/2) 
= (3/2, 5/2, +3/2) 
= (3/2, 5/2, +1/2) 
= (3/2, 5/2, -1/2) 





Magnetic field, B/G 

FIG. 4: (Color online). State-to-state s-wave inelastic cross 
sections for collisions originating in the uppermost hyperfine 
state, /3i8, as a function of magnetic field, for a collision energy 
of 10~^ K (upper panel). The lower panel shows the initial 
and final states energies, color-coded as for the cross sections. 
Solid, dashed and dotted lines represent states with rUs = — 1, 
and +1 respectively. 



with collision energy at fields below ~5 G is simply due 
to a larger probability of tunneling through the outgoing 
centrifugal barriers. 

The behavior of the state-to-state cross sections from 
other states in the = -fl manifold is considerably 
more complicated. Fig. |4] shows the state-to-state cross 
sections for collisions that start in l3ig , which at low field 
is (1/2,3/2,-^3/2) and at high field is [+1, +1, -1/2]. 
Once again the cross sections at low field are mostly gov- 
erned by the kinetic energy release. 

At high field, the strongest transitions are those to 
I3[ = [-1, +1^-1/2] and I3[^ = [0,-^1,-1/2], which are 
driven by the main mechanism with no change in rrii 
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quantum numbers. As before, the transition to = — 1 
is stronger because of the larger kinetic energy release. 
The next strongest are to (i'^ = [—1, +1, +1/2] and j3'^ = 
[0, +1, +1/2], with Amj3 = +1. 



As discussed above, for non-spin-stretched states such 
as /3i8 it is possible to relax nis while conserving m/. 
This is the case for transitions to j3'-j and I3[^, which are 
dominated by L' = and therefore are not centrifugally 
suppressed. However, it is clear from Fig. |4] that the cen- 
trifugally unsuppressed channels are not the dominant 
ones, even at very low field: the hyperfine splittings re- 
lease enough kinetic energy that the main mechanism 
dominates over centrifugally unsuppressed transitions at 
all values of B. 



b. Total s-wave inelastic cross sections. The behav- 
ior of the total s-wave inelastic cross sections with mag- 
netic field, for the three initial states studied above, is 
shown in Fig. [5] for a range of collision energies. If hy- 
perfine interactions are neglected, the quantity cr^^, = 



by 



J, of Eq. (17) is independent of B, AEj3/3> is given 
gslJ-BBAms, and three main regimes are observed 



|24j : (1) at low enough fields, the inelastic cross sections 
flatten out to a zero-field value proportional to E'^; (2) 
as the field increases, cr(^(^ enters a region of B^/^ de- 
pendence, given by the increasing probability of tunnel- 
ing through the c?-wave centrifugal barrier in the domi- 
nant outgoing channel(s); and (3) at high enough fields 
(above about 100 G) the c?-wave centrifugal barriers are 
overcome and o'f^^ again approaches a field-independent 
value, this time proportional to 




10' 10^ 
Magnetic field, B/G 



FIG. 5: (Color online). Total s-wave inelastic cross sections as 
a function of magnetic field for a variety of collision energies 
of 10""^ K (top panel) and 10~^ K (center panel). The states 
are color-coded as in Fig. [i] The dotted line shows the B^^^ 
behavior in regime (2) (see text). 



The inclusion of hyperfine terms modifies both the 
qualitative behavior in regime (1) and the boundaries of 
regime (2). First, at very low fields and collision energies, 
the state-to-state cross sections become nearly constant 
at a field-free value that is much larger than when hy- 
perfine coupling is neglected. The kinetic energy release 
in this region is dominated by Ai?^^', which for some 
outgoing channels does not approach zero as the field 
decreases. The field-free cross section is proportional to 
£■-1/2 at the lowest energies, though deviations from this 
occur at energies high enough that the outgoing energy is 
no longer dominated by AEppi. Increasing the magnetic 
field alters AEppi, particularly above ~10 G, but also 
changes ctJj^, because the spin character of the monomer 
eigenf unctions changes. This leads to a non-power-law 
increase in the cross sections up to the onset of regime 
(2). It may again be noted that, although L' = is pos- 
sible for initial states other than the spin-stretched state, 
the hyperfine splitting at zero field provides enough ki- 
netic energy release for the main mechanism to dominate 
spin relaxation, even though it is suppressed by d-wave 
outgoing barriers. 



4- Prospects for sympathetic cooling 

Trap losses in a static trap are fundamentally caused by 
four phenomena: spin relaxation, background gas colli- 
sions, blackbody radiation and non-adiabatic transitions 
to untrappable states. Non-adiabatic transitions, which 
are one-body transitions that can occur at points in the 
trap where different states are near-degenerate, have im- 
portant consequences for trap design. In particular, it is 
well known for atomic systems that substantial losses can 
occur at the center of magnetic quadrupole traps, where 
the magnetic field is zero [IS]. Near this point, states 
with different values of my are degenerate and the trap- 
ping field varies very fast with position, so that atoms 
can undergo non-adiabatic transitions (Majorana flops 
|49| ) when they pass close to the trap center. For atoms, 
Majorana transitions can be effectively suppressed by ap- 
plying a small bias field (1 G or less) that removes the 
zero-field point. Similar effects have been observed for 
molecules in electrostatic traps |i50i, and can again be 
suppressed by applying a bias field. 

For molecules there are crossings that occur at non- 
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FIG. 6: (Color online) Contour plots of the ratio of elastic to total inelastic cross section as a function of collision energy and 
magnetic field. The panels correspond to calculations excluding (top left) and including (rest) hyperfine terms. The lines show 
the conditions sampled by 99.9% of molecules trapped in the nia = +1 state, in an unbiased trap (solid) and a trap with a bias 
field of B = 1 G (dashed). 



zero field, as shown in Fig. [T] Magnetically trapped NH 
molecules in states /3i3, /3i4 and /3i5 might conceivably 
undergo transitions to untrapped states /3io, /3ii and /3i2 
in the vicinity of crossings that occur between 15 and 
30 G. However, away from the trap center the molecules 
experience a field that varies only slowly as they move, 
and under such circumstances the transition probabilities 
will be very low. We therefore expect that a small bias 
field of 1 G or less will be sufficient to suppress one-body 
losses for NH and other similar molecules. Even if such 
losses do prove significant, the states /3i6, Pn and /3i8 are 
immune to them except near a zero- field point. 

The major loss mechanism in sympathetic cooling 
thus arises from inelastic collisions. Provided that the 
absolute values of the elastic cross sections are large 
enough to provide cooling before the molecules are lost 
to black-body radiation or non-adiabatic transitions, the 
key quantity is the ratio of elastic to total inelastic cross 
sections, which must be greater than about 100 for sym- 
pathetic coohng to proceed. 

Both Mg-NH and NH-NH [FTl ^ collisions may cause 



transitions to untrapped states or release enough energy 
to eject the molecules from the trap. Fig.[6]shows contour 
plots of the ratio of elastic to total inelastic cross section 
as a function of E and B for Mg-NH. The top-left panel 
shows the results when hyperfine terms are neglected, 
and the remainder show the results when hyperfine terms 
are included, for the spin-stretched and two highest-lying 
hyperfine states. 

Trapped NH molecules in state /3 at temperature T 
will be distributed according to a Boltzmann distribution 
with density p given by 



p/po = exp 



kT,T 



(18) 



At any given temperature on the energy axis of Fig. |6j 
only about 0.1% of molecules will experience fields cor- 
responding to energies greater than 6kBT. The diagonal 
lines in Fig. [6] show the maximum fields sampled by over 
99.9% of the molecules trapped in one of the hyperfine 
levels of the rris = +1 manifold. These correspond, re- 
spectively, to a unbiased trap with zero magnetic field 
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at the center (solid lines), and a trap with a bias field 
of 1 G, to prevent Majorana transitions (dashed lines). 
Pre-cooled molecules would enter the trap from the right- 
hand side of the panels in Fig.|6]and then move to the left 
as they are cooled, remaining below the line appropriate 
for the trap in use. 

The ratio of elastic to inelastic cross sections exceeds 
100 at temperatures below 10 niK for all three hyperfine 
states, and is thus favorable for sympathetic cooling to 
work, provided the molecules can be precooled to this 
temperature. This agrees with the hyperfine-free results 
of Wallis and Hutson [53]. At lower temperatures, the 
ratios of elastic to inelastic collisions are not as favor- 
able as in hyperfine-free calculations, but are neverthe- 
less adequate to reach temperatures below 1 /iK. Trap- 
ping molecules in the fully spin-stretched state may be 
particularly advantageous, especially for molecules with 
stronger hyperfine interactions than NH. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the effect of hyperfine interac- 
tions on spin-relaxation collisions of NH with Mg in the 
cold and ultracold regimes. We find that hyperfine in- 
teractions make substantial changes to inelastic collision 
rates at temperatures below about 10 mK and magnetic 
fields below about 20 G. The major effect arises because 
hyperfine interactions modify the kinetic energy released 
in spin-relaxation collisions. When hyperfine interactions 
are neglected, the kinetic energy decreases to zero as the 
field is decreased, but when hyperfine interactions are in- 
cluded the kinetic energy release is significant for most 
transitions even at zero field. For s-wave collisions, the 
kinetic energy release helps overcome the d-wave centrifu- 
gal barriers that suppress spin-relaxation collisions and 
thus leads to larger inelastic cross sections. 

Hyperfine interactions also introduce new mechanisms 
for spin-relaxation collisions. For initial states that are 
not spin-stretched, the cross sections for some of these are 
centrifugally unsuppressed. However, for Mg-NH, where 
the hyperfine interactions are quite weak, the centrifu- 
gally unsuppressed mechanisms make only a small con- 
tribution to total inelastic cross sections at the collision 
energies and fields studied here. It is nevertheless possi- 
ble that centrifugally unsuppressed channels may be im- 
portant in other systems, with either stronger hyperfine 
interactions or weaker competing spin relaxation mecha- 



nisms. 

The most important hyperfine effects for Mg-NH arise 
from the Fermi contact interactions. These determine 
both the composition of the n = states in terms of 
spin functions and the low-field energy level pattern (and 
hence the kinetic energy release). Other hyperfine terms 
have only very small effects for transitions between n = 
states. Indeed, we have repeated the calculations of 
the state-to-state cross sections including only the Fermi 
contact interactions and obtain almost identical results. 

Our results for Mg-NH('^I]^) may be compared with 
those of Tscherbul et al. [53] for He-YbF(2E). For YbF 
the main mechanism of electron spin relaxation consid- 
ered here, driven by the electron spin-spin coupling, does 
not exist and is replaced by a higher-order and much 
weaker mechanism driven by the spin-rotation interac- 
tion. Under these circumstances, combined electron and 
nuclear spin relaxation, driven by the electron-nuclear 
dipolar interaction, is in relative terms much more im- 
portant. However, Tscherbul et al. did not focus on the 
regime where hyperfine energies make important contri- 
butions to the kinetic energy release. 

We have considered the prospects for sympathetic cool- 
ing of NH by Mg, which were previously explored in 
hyperfine-free calculations by Wallis and Hutson [24] . We 
have calculated the ratio of the elastic to inelastic cross 
section as a function of energy and magnetic fields for sev- 
eral magnetically trappable hyperfine states of NH. Even 
though hyperfine interactions increase inelastic cross sec- 
tions at low energies and magnetic fields, the ratio re- 
mains high enough for sympathetic cooling to proceed if 
the NH molecules can be pre-cooled to about 10 mK. 

Molecular hyperfine interactions are also likely to be 
important in developing techniques for controlling ultra- 
cold molecules. The low-lying excited states afforded 
by hyperfine splittings can support near-threshold lev- 
els that will produce magnetically tunable Feshbach res- 
onances. Once molecules such as NH have been cooled to 
the ultracold regime, it will be possible to use such reso- 
nances both to control collisions by adjusting the scatter- 
ing length and to create polyatomic molecules by mag- 
netoassociation, as has already been achieved for alkali- 
metal atoms [T31 [13] . 
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